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Abstract 

In various extensions of the Standard Model (SM) tree level non-leptonic decays of hadrons 
receive contributions from new heavy gauge bosons and scalars. Prominent examples 
are the right-handed W bosons in left-right symmetric models and charged Higgs {H^) 
particles in models with extended scalar sector like two Higgs doublet models and super- 
symmetric models. Even in the case of decays with four different quark flavours involved, 
to which penguin operators cannot contribute, twenty linearly independent operators, in- 
stead of two in the SM, have to be considered. Anticipating the important role of such 
decays at the LHCb, KEKB and Super-B in Rome and having in mind future improved 
lattice computations, we complete the existing NLO QCD formulae for these processes 
by calculating 0(a s ) corrections to matching conditions for the Wilson coefficients of all 
contributing operators in the NDR-MS scheme. This allows to reduce certain unphysical 
scale and renormalization scheme dependences in the existing NLO calculations. Our 
results can also be applied to models with tree-level heavy neutral gauge boson and scalar 
exchanges in AF = 1 transitions and constitute an important part of NLO analyses of 
those non-leptonic decays to which also penguin operators contribute. 
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1 Introduction 

In the Standard Model (SM) the non-leptonic AF = 1 decays of mesons are governed by 
the (V — A) x (V — A) structure of the leading four-quark operators originating in the 
tree-level exchanges. If all the four flavours of the participating quarks are different 
from each other the only possible diagrams contributing to these decays in the SM and in 
any of its extensions are the current-current ones: penguin diagrams are absent. Decays 
of this type are theoretically cleaner than the ones in which also penguin diagrams and 
penguin operators contribute. As such they are well suited for the determination of the 
CKM parameters, in particular the angles 7 and j3 in the unitarity triangle jTJ |2j [3]. 

While non-leptonic decays are subject to significant non-perturbative uncertainties orig- 
inating in hadronic matrix elements of four-quark operators, the QCD factorization ap- 
proach to non-leptonic two-body decays [I] combined with advanced lattice calculations 
could one day promote non-leptonic two-body meson decays to precise tools in testing 
the SM and its extensions. In these studies renormalization group short distance QCD 
effects play an important role. In the SM they are known including the NLO corrections 
and in a few processes at the NNLO level. An up-to-date review can be found in [5]. 

Beyond the SM new local four-quark operators with different Dirac structures can be 
generated. The simplest example are (V + A) x (V + A) operators originating in the 
exchange of W^' gauge bosons in the left-right symmetric models. However, also right- 
handed (RH) couplings of the SM W gauge bosons can be generated in various extensions 
of the SM like left-right symmetric models and generally also models with vectorial heavy 
quarks that mix with the SM chiral quarks. In this case also (V — A) x (V + A) operators 
contribute. The latter operators generate through QCD corrections (S — F) x (S + F) 
operators present also in models with charged (-ff ± ) Higgs particles. In the latter models 
also (S ± P) x (S ± P) operators are present. Needless to say all these statements also 
apply to neutral gauge bosons and scalars mediating AF = 1 transitions. 

The full set of twenty linearly independent dimension six four-quark operators with four 
different flavours in all extensions of the SM has been listed in [HI [7J |8] , where also two- 
loop QCD anomalous dimensions of these operators have been calculated. However, the 
full NLO QCD renormalization group analysis of non-leptonic decays requires also the 
calculation of 0(a s ) corrections to matching conditions for the Wilson coefficients of the 
operators in question. While such corrections are known within the SM [9j [10] at the NLO 
level, to our knowledge a complete analysis of these corrections including all operators in 
any extension of the SM is absent in the literature. 

In a recent paper [TT] we have calculated 0(a s ) corrections to matching conditions for the 
Wilson coefficients relevant for AF = 2 processes mediated by heavy colourless neutral 
gauge bosons and scalars reducing thereby certain unphysical scale and renormalization 
scheme dependences present in the absence of such corrections. Similar unphysical scale 
and renormalization scheme dependences are present in the absence of 0(a s ) matching 
corrections in AF = 1 amplitudes generated by tree-level gauge boson and scalar ex- 
changes and it is desirable to reduce them as well. 

The main goal of our paper is the calculation of 0(a s ) corrections to matching conditions 
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for the Wilson coefficients of all dimension six four-quark operators with four different 
flavours contributing to AF = 1 decays mediated by colourless gauge bosons and scalars 
in the NDR-MS scheme. As the two-loop anomalous dimensions for these operators have 
been already calculated in this scheme in [8] our calculations complete the NLO QCD 
analysis of the decays in question. 

Our paper is organized as follows. In Section [2] we recall the general structure of the 
effective Hamiltonians for AF = 1 processes in question and we give the full list of four- 
fermion operators that contribute to these transitions. Subsequently we collect their one- 
and two-loop anomalous dimension matrices. In Section [3] we calculate 0(a s ) corrections 
to the amplitudes in the full theory and in Section H] the corresponding results for the 
matrix elements of operators are presented. This allows us in Section [5] to present the 
Wilson coefficients of all twenty operators including 0(a s ) corrections. In Section [6] 
combining our results with the known renormalization group evolution matrices we arrive 
at a complete NLO formulae for the Wilson coefficients of the involved operators. In 
Section [7] we demonstrate the scale independence of the physical amplitudes analytically 
and in Section [8] we investigate the removal of this scale dependence numerically. We 
conclude with a brief summary in Section [91 



2 Theoretical Framework 
2.1 Local Operators 

While in the SM only two current-current operators contribute to each AF = 1 transition, 
the list of current-current operators beyond the SM is much longer. As in j8] we choose 
the operators in such a manner that all the four flavours they contain are different. In 
such a case, the only possible diagrams are the current-current ones. In what follows we 
will fix the four flavours to be b, u, c, d but other choices are clearly possible without 
changing our results. 

Twenty linearly independent operators can be built out of four different quark fields. 
They can be split into eight separate sectors, between which there is no mixing. The 
operators belonging to the first two sectors (VLL, VLR), that are relevant for gauge 
boson contributions, are given as follows 

QVLL = ^p LU ^^YP L d a ), (la) 
Q2 LL = (b a ltJi P L u a )(c^P L d^, (lb) 

Q^ LR = (bVLuP)(c^P R d a ), (lc) 
QVLR = (6« 7M P £W «)(g0yp fl A (Id) 
where a, ft denote quark colours. In the case of scalar contributions the following operators 
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have to be considered: 

gSLR = (b»p LU ^^P R d a ), (2a) 

gSLR = (l»p LU *)(cPP R dP), (2b) 



gSLL = ^p LU P)^P L d a l (3a) 

gSLL = (l»p LU ^^P L d f3 ), (3b) 

gf L = (b a a^P L u^a^P L d a ), (3c) 

Qf L = (6V^P i iO(c / V*P i d /3 ). (3d) 



The operators belonging to the four remaining sectors (VRR, VRL, SRL and SRR) are 
obtained from the above by interchanging Pl and Pr. Obviously, it is sufficient to calcu- 
late the Wilson coefficients only for the VLL, VLR, SLR and SLL sectors. The "mirror" 
operators in the VRR, VRL, SRL and SRR sectors will have exactly the same properties 
under QCD renormalization. 

The two-loop anomalous dimensions for these operators have been calculated in the NDR- 
MS scheme in [8] with a particular choice of the evanescent operators. As discussed 
there, while these operators are essential for the correct evaluation of two-loop matrix 
elements, the virtue of the formulation of the NDR-MS scheme introduced in [10] is 
that the evanescent operators defined in this scheme influence only two-loop anomalous 
dimensions. By definition they do not contribute to the matching and to the finite gluon 
corrections to the matrix elements of renormalized operators calculated by us. They are 
simply subtracted away in the process of renormalization. This issue is summarized in 
Section 6.9.4 of [12], where further references can be found. A very important paper in 
this respect is also the one of Herrlich and Nierste [13], where this issue is discussed in 
full generality. Therefore effectively the one-loop calculations presented here and in [UJ 
are based on the projections of various Dirac structures on physical operators that are 
consistent with [8j [10] but otherwise the evanescent operators can be dropped from the 
beginning. 



2.2 Renormalization Group Functions 

2.2.1 Running of QCD Coupling and Running Masses 

For the complete NLO renormalization group analysis we need a number of renormaliza- 
tion group functions that we recall here for completeness. 

In particular we need the QCD (3 function at the two-loop level 
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where 



UiV-2/ 34 ^,10 C F = ^1 

H 3 3 3 2N 



(5) 



/ is number of flavours and N the number of colours. 

Similarly, the two-loop expression for the quark mass anomalous dimension can be written 

2 



as 



A 1 ) 



a. 



where 



7i 0) = QC F t£) = c F (sc F + ^iV - H/) . 



(6) 
(7) 



2.2.2 One-Loop and Two-Loop Anomalous Dimension Matrices of Operators 

The most important ingredients of any renormalization group analysis in weak decays are 
the anomalous dimension matrices that we define in general form as follows 



l(oc s 



1 An 1 \AnJ ' 



(8) 



In particular the one-loop anomalous dimension matrices 7^ will play important role in 
our discussion of the removal of the unphysical fi dependences at NLO in Section 

The one-loop matrices for all operators considered in this paper read 



7 



(O)VLL 



-4 6 

N u 

6 

u N 



(9) 



7 



(O)VLR 



-6N+j? 



-6 



N 



(10) 



^(O)SLR _ / N 6 

1 -6AT + |p 



(11) 



7 



(O)SLL 



V 



_6_ 

N 




f + 247V 



-6iV + f 



24 

_48 

N 



N_ 1_ 

2 AT 
1 

-1. - 4iV 

AT 4:1 y 



1 

2 

_ J_ 

AT 

6 



\ 



(12) 



48 -f 2iV-f/ 

The two-loop matrices for all operators in the NDR-MS scheme used in this paper read 
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Lsj 



_ 22 _ 57 __2_f _ 19 /U _|_ 39 _i_ 2 x 

~(1)VLL _ / 3 2iV 2 3Af^ 6 iV JV 3^ \ /i o\ 

> _ i£ /V _i_ 39 i 2 ^ _ 22 _ _57_ 2_ x J > ( ° 1 

6 JV 3^ 3 2JV 2 3JV 

203 7\r2 , 479 , Jj3_ , 10 at x 22 f 

6 1 V " r 6 ^ 2W 2 ^ 3 JV J 37V J 

_ 100 /V _i_ _3_ _|_ 22 L^n 15 _ _22_ x 

3 JV 3 ^ 6" 1 " 2JV 2 3JV-/ ■ 

137 I 15 _ _22_ x _ 10U at I 2 22 f 

(l)SLR _ / 6 ~r 2JV 2 3JV-> 3 TV 3 » 



19 7vr _|_ 39 _|_ 2 J? 22 57 2_ x 

~ 6 JV 3 1 * 3 27V 2 3JV^ 

_ 203 at2 , 479 , JJ_ i 10 at f _ 22 x _ 71 at _ 18 , a f 

-(l)VLR / 6 ly 6 ~ r 2N 2 ^3 ly J 3N J 2 ly N ' \ (-ia\ 

"1 100 at I 3 I 22 x 137 ,15 22 x J ' 1 ! 1 1 



" \ _71at_ 18 +4f 203 T\j2 _i_ 479 j 15 i_ 10 at x 22_ x J ' 

2 iv n ^ 6 iV ^ 6 ^ 2JV 2 ^ 3 iV ^ 3/V^/ 

(l)SLL _ _ JV 2 , 148 _ _L07_ QATf 10 f 
'11 2 "T 3 2/V 2 ZJV ^ 3JV->> 

(l)SLL _ _178 m | 64 , 16 f 

712 - 3 JV /V ~+~ 3 J ' 

(l)SLL _ 107/V 2 _Zl_J__J_/Vf+X 

'13 - 36 iV 18 /V 2 18 iV ^ + 9iV' 

(l)SLL = _ 109 at , A _ X 
'14 36 JV 18' 

7 « SLL = -26/V + f , 

(l)SLL _ _203 at2 , 28 107 , 10 at X 10 r 

722 - — iv +T 2F f I JV J 3JV-/' 

(l)SLL _ 89 at , 2_ _ 1 x 

723 - 18 JV + JV 9/' 

(l)SLL _ _53 _ _4_ , J_x 

'24 is JV 2 9JV ^ ' 

(l)SLL _ 676 at2 _ 1880 _ 320 _ 88 a t f , 176 x 

731 - 3 JV 3 JV 2 3 + 3JV-/' 

(l)SLL _ 820 at I 448 _ 88 f 

732 3 JV 3 •> ' 

(l)SLL _ _2^ 7V 2_116,^l_,22 A T f ,^_ f 

'33 - 18 JV 9 2JV 2 9 ^ SN^i 

(l)SLL _ 50 at 8 x 

'34 - 3 JV 3-7> 

488 at _i_ 416 176 x 

3 JV 3 •> ' 



7 (1,SLL 



(l)SLL _ 776 320 , 176 x 
'42 3 JV 2 " 3JV J ' 



(l)SLL _ 22 at 40 _j_ 8 £ 

43 - 3 JV JV 3-/' 

343 A T2 + 28 + ^l_26 A T f+ ^_ f QgN 
18 9 2JV 2 9 1S J ^ 9NJ- 



7if LL 
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2.3 Effective Hamiltonian 

The effective Hamiltonian for AF = 1 transitions can be written in a general form as 
followa3 

n A e t l ^ K^CMQi , (17) 

i,a 

where Q\ are the operators given in Eqs. ([I]), and ([3]) and Cf (/i) their Wilson coeffi- 
cients evaluated at a scale \i at which the hadronic matrix elements are evaluated. The 
scale /i can be the low energy scale /i^ at which actual lattice calculations are performed 
or any other scale, in particular the matching scale /j, in . In this case the matrix elements 
are obtained by evolving lattice results by means of renormalization group (RG) equa- 
tions from /i£ to /ii n . The resulting matrix elements of the effective Hamiltonian, that are 
directly related to decay amplitudes, can then be written generally as follows: 

CH A c i =1 } = nJ2Q(^)(Qt(^)} ■ (18) 

i,a 

The overall factor k in our analysis depends on the exchanged boson (vector or scalar) 
and will be chosen such that for non- vanishing Wilson coefficients Cf(ni Q ) = 1 in LO. 

Evidently the matrix elements (QKfim)) depend on the matching scale /i in and also on the 
renormalization scheme for operators. In order to remove these unphysical dependences 
from physical amplitudes one has to calculate 0(a s ) corrections to Cf(//j n ). This is the 
goal of Sections [3HSJ 

2.4 Procedure for Matching 

Before entering the details let us recall that the calculations of 0(a s ) QCD corrections 
to Wilson coefficients are by now standard and have been described in several papers. 
In particular in [T3| [T2] in the case of the SM. In [11] this procedure has been used to 
calculate 0(a s ) corrections to the Wilson coefficients of operators contributing to AF = 2 
FCNC processes mediated by tree level neutral gauge boson and scalar exchanges. Here 
we recall briefly this procedure that we will use in the case of AF = 1 decays. 

Step 1 

We first calculated the amplitudes in the full theory. This amounts to the calculation of 
the diagrams in Figs. [I]and[2j in the case of a gauge boson exchange and a scalar exchange, 
respectively. In the presence of massless gluons one encounters infrared divergences. We 
have regulated these divergences by a common external momentum p with — p 2 > for all 
external massless fields as we did in [TTJ. The ultraviolet divergences present in diagrams 
in Fig. [T] and |2] with gluon corrections to vertices have been regulated using dimensional 
regularization with anti-commuting 75 in 4 — 2e dimensions. 



1 In what follows we drop for simplicity h.c. 
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Figure 1: Tree level diagram and one-loop QCD corrections to AB = 1 transition 
mediated by a gauge boson in the full theory. Mirror diagrams are not shown. 




Figure 2: Tree level diagram and one-loop QCD corrections to AB = 1 transition mediated 
by a scalar particle in the full theory. Mirror diagrams are not shown. 



Step 2 

We have calculated the matrix elements of contributing operators by evaluating the di- 
agrams in Fig. [3] making the same assumptions about the external fields as in the first 
step. In contrast to step 1 one has to renormalize the operators. This we do in the MS 
renormalization scheme with anti-commuting 75, which corresponds to the NDR scheme 
used also in [8] and [15] • 
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Figure 3: Leading order and one-loop diagrams in the effective theory. Mirror diagrams 
are not shown. 



Step 3 

We finally inserted the results of the two steps above into a formula like the one in Eq. ( fl8i) 
and comparing the coefficients of operators appearing on the l.h.s (full theory) and r.h.s 
(effective theory) we found the coefficients C?(/ii n ). As these coefficients cannot depend 
on the infrared behaviour of the theory, the dependences on p 2 found in the first two steps 
have to cancel each other in the evaluation of C^(fii n ). Indeed we verified this explicitly 
The interested reader can do this as well by inspecting our intermediate results that we 
present in Sections |3] and HI respectively. 

Very often in analyses of new physics contributions the overall factor in front of the sum 
in Eq. ( !T8|) is chosen as in the SM. However, in our analysis it will be more convenient to 
use in each case the normalization in which the Wilson coefficient of the leading operator 
evaluated at the matching scale is equal to unity in the absence of QCD corrections. In 
this manner the applications of our formulae in various new physics (NP) models will be 
facilitated. 

After this preparation we are ready to present our calculations in three steps in question. 

3 Amplitudes in the Full Theory 

Using the Feynman rules in Fig. H]we find for the colourless gauge boson exchange (after 
quark wave function renormalization) 
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i^Saf, [A l [(A)P L + A%(A)P R ] 
i5 a ? [A*l(H)P L + A%(H)P R ] 



Figure 4: Feynman rules for colourless charged gauge boson A + with mass Ma, and 
charged colourless scalar particle H + with mass Mh, where i (j) denotes an up-type (down- 
type) quark flavour with charge +| (—\) and a, (3 are colour indices. 



A 



VLL 



(At(A)Y Af(A) 



M\ 



O.s 



1 + 2C F -^ lo, 

47T 



/r 



-p 



-oVLL 
2 / ^2 



a s fl . M 2 A 
+ 47U + 1 ° g ^ 



(19) 



Q v 2 hh - 3QI 



N 



VLL 



A 



VLR 



Mi 



1 + 2C F -^ lo, 



Air 



-p 



n VLR 
2 / ^2 



+ 



4vr \2 



+ log 



_p2 



(20) 



gy LR + 3g x v 



AT 



VLR 



For the colourless scalar exchange (after quark wave function and quark mass renormal- 
izations) we find 



SLR 



(Af(H)y A%(H) 
Ml 



a 



1 + 8C F -^ ( 1 + log 
Air 



—p 2 



Ql LR 



(21) 



A! 



SLL 



(At(H))*A?{H) 



2M 2 H 



a 



1 + 8C F -^ 1 + lo 



+ 



Att V 2 



+ log 



4-7T 



— 



1 

2iV 



n SLL ^ n SLL 
V4 



(22) 



4 Matrix Elements of Operators 

Calculating the diagrams in Fig. [3] we find after quark wave function renormalization and 
operator renormalization 
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(Qi ) 



<gy LL > 



VLR\ 



{Ql 



<gy LR > 



(Qf R ) = 



(Ql LL ) 



(Ql LL ) 



(Qf L ) 



-p 2 

,,2 1 



-gy LR 

TV 



3 n VLL 

iV^ 2 
3Q^ LR 



- 3gy LL 

■ 3Q^ LL 



i + 2 ^(i + io g ^)]gy LR 



(Qf R ) 



(Qf L ) = 



(23) 
(24) 
(25) 

(26) 



— » 2 



1 + 8C F ^ (l + log-^ 

47T \ — p 2 



; + log — 



Q2 + 47 



-p 2 

3 n VLR 

TV 



3 n SLR 

- 3g^ LR 



+ 3g 2 SLR 



1 + 2C F ^ 

47T 



47T \b —% 



-p L 



--Qf LL 
7V Vl 



+ -± 2 + logA 
Air \ —p- 



,2 x 

^2 



'9 _ /V 2 1 

- — — g? LL - -gf L 

4AT V3 4 V4 



47T \ —J. 



2 * 
|2 



-p 2 

1 



gf L 



47T \ ^ — 

1 + 47 3iV ( l0g ^ + ^ 



^QfL _ IgSLL 



(27) 
(28) 

+ 3gf L 

(29) 
(30) 



+ 4 N 



a s _ / 2 

+ 



47T 



■12 



6iV 2 



«3 ll -S 3 (} +io ^K l 



(31) 



/oSLL 
2 / ^2 
'V J 



, qi- + 24^(i + log .£) [±g§"--«r 



47T 



2 n SLL 

TV 



- 2Qf L l 



(32) 
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ff 



We remark that for the matching performed in this paper only QCD corrections to the 
matrix elements (Qy LL ), (Qy LR ), (Q| LR ) an d (Ql LL ) are required. The QCD corrections 
to the matrix elements of the remaining operators would enter the NLO analysis if the 
exchanged gauge bosons and scalars were coloured. They would also enter our analy- 
sis if we included 0(a 2 ) corrections. We include these additional matrix elements for 
completeness. 



5 Results for the Wilson Coefficients 

In what follows we will list the general structure of the effective Hamiltonian in each case 
and subsequently we will list our results for the Wilson coefficients (we drop h.c. in what 
follows). 



5.1 Colourless gauge boson 



{Af(A)Y Al d (A) 



?4T se = LK iL [C 1 VLL (//)gy LL + C 2 VLL (/i)Q 2 VLL ] 



M A 



4vr V /i 2 2 

4tt V N ° g fx 2 2N J " ' 4tt V ^ V 2 6 



(33) 



+ ( A L(A))^(A)) 2 ^ VLR( ^ )gVLR + ^vLR^gyLRj + L ^ R 

We find for an arbitrary number of colours N 

CrM ^(_ 31og M + ^, (34) 

vi t v a s /3 , Ml 11 \ a s / Ml 11\ . . 

^' LR M = ? i f31og^ + iV (36) 



CS™M = 1 + ^ f 4 ■<* # - i V 1 + £ (~ log # - i ] ■ (37) 
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5.2 Colourless scalar 



- (Af( ^ A ' (A) [Ci LL MQT L + Cf^)Qf L ] (38) 

- (Af(g ^ 2 A ^ (iJ) [^ SLR (/i)Qf R + C^)Qf R ] + L o i? 
We find for an arbitrary number of colours N 

Cf^fji) = 3^ , (39) 

(40) 

Cf i (A*) = 0, (41) 
C 2 SLL (/i) = l, (42) 

^•)-£(^9-5r)-£(N9-s)- 

We emphasize that the Wilson coefficients C" of the "mirror" operators (Pl -h- Pr) as 
defined by us are equal to the ones presented above. The formulae presented in this 
section are the main results of our paper. 



6 Master Formulae for NLO Wilson Coefficient func- 
tions 

With these results at hand and the known one-loop and two-loop anomalous dimension 
matrices that we have listed in Section [2] we can complete the NLO renormalization 
group analysis which can give us the Wilson coefficients at the low energy scales at which 
hadronic matrix elements are evaluated. In the case of the B decays the final result for 
the Wilson coefficients in each of the four sectors of operators with no mixing between 
different sectors can be written as follows: 



C(fl b ) = U(fl b , /X in )C(/Xin), 



(45) 
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where U are evolution matrices and C column vectors. For each VLL (VRR), VLR (VRL), 
SLR (SRL) sector U is a 2 x 2 matrix, while it is a 4 x 4 matrix in the SLL (SRR) sector. 
Similarly for each VLL (VRR), VLR (VRL), SLR (SRL) sector C is a two-dimensional 
column vector, while it is a four-dimensional one in the SLL (SRR) sector. 

General formulae for the evolution matrix at the NLO level have been derived in [16]. In 
order to make our paper self-contained we recall these formulae in what follows: 

tWta) = (l + ^-j) ^(Min) (l " ^-j) ■ (46) 

Here £A°) is the evolution matrix in leading logarithmic approximation and the matrix J 
expresses the NLO corrections to this evolution. We have 



C/ (0) (/^in)=t> 

where V diagonalizes 



2d () 



V- 1 , (47) 

D 



^g) = y-i^ )Ty (48) 

and 7^°- ) is the vector containing the diagonal elements of the diagonal matrix 7^. 
If we define 

G = V-i^O-nv (49) 

and a matrix H whose elements are 



(0) Pi Gij 

llij — Uif]' 

the matrix J is given by 



B« = - ——^ ^ , (50) 



J = VHV~ 1 . (51) 

We next write our results for the Wilson coefficients at the matching scale in a general 
form as follows 

d(lM*) =C - ^p^C, . (52) 
Finally combining these initial values with the evolution matrix f l4"6]) we obtain 

C(n b ) = (l + ^-Jj U {0 \n b ,n- m ) (l - (C 1 + J(7 )) • (53) 



We recall that when using this evolution down to low energy scales one has to insert the 
correct number of effective flavours. As this procedure is by now standard, we refer to 
Section HIE in [17] for details. 

We end this section by recalling certain features of the fundamental formula (|53l) : 
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• The renormalization scheme dependence of the matrix J on the left-hand side of 
the LO evolution matrix is cancelled by the one of hadronic matrix elements. 

• This scheme dependence on the right-hand side of the LO evolution matrix is can- 
celled by the one of C\ calculated by us: C\ + JCq is renormalization scheme inde- 
pendent. 

• The dependence on the precise choice of the scale fib in the evolution matrix is 
cancelled by the one present in the hadronic matrix elements. In the case of AF = 2 
transitions in the SM, in which only a single operator is present this allows to 
introduce renormalization scheme and renormalization scale invariant parameters 
like B K . 

• The dependence on the precise choice of the scale fi- in in the evolution matrix is 
cancelled by the logarithmic terms in C\ that we calculated in this paper. 

We will now look in more details at the last issue. 



7 Renormalization Scale Dependence 

One of the main virtues of our calculation of 0(a s ) corrections to Wilson coefficients 
at the high energy matching scale fi- m is the cancellation of the fi- in dependence of the 
renormalization group evolution matrix by the fi- m dependence of the Wilson coefficients in 
question. This cancellation requires particular values of the coefficients of the \og(M 2 /fif n ) 
in Cj(/ij n ) where M stands for the mass of a heavy gauge boson or heavy scalar involved. 
As this cancellation constitutes an important test of our results it is useful to derive a 
general condition on the coefficients of log(M 2 / ft 2 n ) in Cifa n ). 

To this end let us look at the evolution matrix in Eq. (1431) . Expanding this matrix around 
the two fixed scales m& and M keeping only the logarithmic terms one obtains 

Ufa, A*0 = (1 + log U(m b , M) [1 + — log j , (54) 

where 7 (0) is the coefficient of a s in the one loop anomalous dimension matrix that de- 
scribes the mixing of operators. The 7 (0) matrices for VLL (VRR), VLR (VRL), SLR 
(SRL) and SLL (SRR) sectors have been collected in Section O Note that it is ^°^ T and 
not y(°> that enters Eq. ( l5"4"j) . Moreover, in the study of the fi m dependence in the case of 
the scalar exchange one has to take into account that in this case the m 2 (fi in ) dependence 
is hidden in the coefficients (A^ //? (#))*A^ //? (#). 

Considering then the cases of colourless gauge bosons and scalars we find that the following 
quantities should be fi- m - independent: 

i? gauge = Ufa fi in )Cfa n ) , (55a) 



if-lar = ^Cfa^m'fa^ 



(55b) 
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C(^)=C -^£^Klog 1 -^ r , (56) 



For each VLL (VRR), VLR (VRL), SLR (SRL) sector C(/z) is a two-dimensional column 
vector, while it is four-dimensional for the SLL (SRR) sector. 

We write next in each case 

4-7T 

where we suppressed ^ m independent 0(a s ) terms and 

mVin) = m\M) (l + ^^ 7 i 0) log ^) , (57) 
with governing the scale dependence of quark masses in QCD. 

Imposing Eq. (j33|) . the conditions for K to ensure /x in independence of resulting amplitudes 
in these two cases read 

^ gauge = 1—Co , (58a) 

K sc ^ = 1^-00 + ^00. (58b) 

Thus the coefficients of logarithms in (7(/Xi n ) can be found without the calculation of 
the loop diagrams in Fig. [T]|3] but formulae in Eq. (1351) serve as a useful check of our 
results for logarithmic terms. These terms are renormalization scheme independent and 
while cancelling the /i; n dependence of {/(//&, fMn) in perturbation theory cannot remove 
its renormalization scheme dependence at the NLO level. To this end as discussed in 
Section [BJ the 0(a s ) non- logarithmic terms have to be calculated which constitutes the 
main new result of our paper. 

Inserting the formulae for the one- loop anomalous dimension matrices and 7m , that we 
listed in Section [21 into Eq. (135]) one can verify that the resulting coefficients K sange and 
^-scalar e q U£ j ^he coefficients of the logarithmic terms calculated by us. This implies that 
the inclusion of 0{pc s ) corrections in question remove the unphysical dependence on the 
precise value of the matching scale. 

The manner in which the /ii n -dependence is removed resembles the one which we encoun- 
tered in the case of AF = 2 transitions HTl: 



In the gauge boson case the /^-dependence of U /ii n ) can only be cancelled by 
the corrections calculated by us. 

The case of scalar exchange with SLR couplings is quite different. Here the /in- 
dependence of U(fJLb, /iin) is totally cancelled by the one of the m 2 (/i; n ) so that even 
without our corrections the amplitudes £1X6 /J'in independent. The role of our calcu- 
lation in this case is then the removal of the renormalization scheme dependence. 

In the case of SLL operators both the running quark masses and the corrections 
calculated by us are required for the removal of the unphysical matching scale de- 
pendence present in LO. 



8 Numerical Analysis 



16 



8 Numerical Analysis 

We will now compute the size of unphysical //^-dependence present in the LO expressions 
and we will demonstrate their reduction after the inclusion of 0(a s ) corrections. 

Compared to our analysis of AF = 2 processes in [11] a numerical analysis of left-over 
unphysical scale dependences in the decay amplitudes is complicated by three facts: 

• The analogues of Pf factors [15] that summarize the renormalization group effects 
between high energy and low energy scales and include also the values of hadronic 
matrix elements do not exist in the literature for decays discussed here although 
obviously complete NLO analyses for the VLL case are known. 

• The hadronic matrix elements of new operators contain larger theoretical uncertain- 
ties than in the case of AF = 2 decays. 

• The increased number of operators relative to the ones in the case of AF = 2 
transitions makes the analysis in questions very involved. 

Leaving a detailed analysis in a concrete model for the future, we nevertheless illustrate 
the size of unphysical /z in -dependence present in the LO expressions in the sectors VLL, 
VLR and SLL. In the SLR sector our corrections have no impact on the cancellation of 
the /ij n dependence as we discussed above. 

We will solve the problem of poorly known hadronic elements in the following manner. 
We will go to the operator basis in which the one-loop anomalous dimension matrices are 
diagonal. At LO in each of the VLL, VLR and SLR sectors there are then two-operators 
that do not mix with each other and with other operators of other sectors. There are four 
such operators in the SLL sector. This property remains true at NLO only in the VLL 
sector as only in this sector the one-loop and two-loop anomalous dimension matrices are 
diagonalized by the same matrix V. For VLR and also scalar sectors one has to use the 
general formulae presented in Section |6j 

Yet, as we have seen in the previous section in order to study the cancellation of the 
/independence present in the LO expressions it is sufficient to study the evolution matrix 
at the leading order and keep only the leading logarithms in Diagonalizing then the 
matrices in all four sectors we can factor out the hadronic matrix element in each 
case so that the issue of the study of the /ti n -dependence as expected can be investigated 
transparently without any hadronic uncertainties. 

Denoting by C± the Wilson coefficients of the operators Q± corresponding to the eigen- 
values of the two anomalous dimension matrices a = VLL, VLR and normalizing C± so 
that they are equal unity at the matching scale at LO we find: 



VLL Case 



Q 



,VLL 




i VLL 



'VLL _ ^-yVLL i ^VLL 
+ — Lj 2 T 1 , 




(59) 



2 
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Q- 



VLL 



gyLL 



Ql 



VLL 



c 



VLL 



VLR Case 



£<VLL 



VLL 



(0)- 
7vLL 



6 fi 

•iv" 6 - 



(60) 







VLR 
+ 







VLR 
2 



VLR 



N 



C 



VLR 



G 



VLR 



(0)+ 
7VLR 



N' 



(61) 



VLR 



VLR 



iV 



c 



VLR 



C-VLR + NC ™, 



(0)- 
7VLR 



N 



-6N. 



(62) 



Here 7i°^ ± denote the anomalous dimensions of the operators Q±. Then the quantities 
(a = VLL, VLR) 



R. 



± 



C^s (/^in 

« S (M)J 



2/3 



1 + 



47T 



M 2 



(63) 



where denote the non-logarithmic 0(a s ) corrections in the new basis should be 
independent at 0(a s ). Here M is the mass of the exchanged gauge boson that we will set 
to 1 TeV in our numerical calculations. 

We should emphasize that in addition to 0(a s ) corrections proportional to there are 
also corrections at this order that come from the evolution matrix. We have denoted them 
by J in Section [HI We do not include them in our analysis in order to exhibit the size 
of the corrections calculated here but they have to be included in any phenomenological 
analysis in order to obtain renormalization scheme independent results. 

From our discussion of the previous section and the diagonalization performed here we 
find that 



la 



a = VLL, VLR 



(64) 



11(N 



VLL 



VLR 



2N 



11 

y ' 



2N 



1 

6 



11(JV + 1) 



VLL 



VLR 



2N 



N 2 -l 



22 
T 



2N 



4 
3 



(65a) 



(65b) 



and indeed then the should be equal unity after the inclusion of only the 0(a s ) 
logarithmic corrections calculated here. The departure from unity at NLO signals 
the presence of terms. 

In Fig. [S]we plot R^ for a = VLL, VLR as functions of the matching scales setting as an 
example the masses of gauge bosons to 1 TeV for three cases: Without the 0(a s ) cor- 
rections in the Wilson coefficients (dashed blue line), with the contribution proportional 
to K± (dotted green line) and including both logarithmic and non- logarithmic 
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500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 



Pin Pin 

Figure 5: The quantities Rf (a = VLL, VLR) defined in Eq. $63) as a function of \ii n for 
M = 1 TeV. The LO result (removing contribution proportional to and r^ ) is shown 
by the dashed blue line, the dotted green line is the NLO result including only logarithmic 
0(a s ) corrections and the solid red line shows the NLO result including both and 



corrections (solid red line). While the dashed blue lines exhibit a significant /i in depen- 
dence, the dotted green lines and solid red lines stay nearly constant over the considered 
/ii n range. The dotted green lines that include only logarithmic 0(a s ) corrections are 
equal to 1 at fi in = 1 TeV as expected. For the solid red lines a small shift relative to the 
dotted green ones occurs which is due to the non-logarithmic corrections r^. Only in the 
VLL sector which has larger than the VLR sector a slight n\ n dependence occurs. This 
dependence can only be cancelled by NNLO corrections. 

SLL Case 

The diagonalization in the SLL case is a bit more involved since four operators are present. 
Furthermore, in the scalar case we also have to take into account the /ij n dependence of 
the quark masses (see Eq. ( 155]) ). Thus instead of Eq. ( 163|) the following quantities 



T33 

SLL 



a s (M) 



J0)j 
SLL 

2/3 



47T 



"^SLL 



+ ^SLL 



ra 2 (/x in ) 
m 2 (M) ' 



±±,± 
(66) 



should be /i in independent at 0(a s ). For simplicity we set in the following N = 3 in order 
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Figure 6: The quantities R 3 SLL (j = ±±, ±) as a function of [ii n for M = 1 TeV. The LO 
result (removing contribution proportional to K 3 SLL and r 3 SLL ) is shown by the dashed blue 
line, the dotted green line is the NLO result including only logarithmic 0(a s ) corrections 



K J SLL and the solid red line shows the NLO result including both K J SLL and r{ 



SLL ' 



to shorten the formulae. We find 

(o)i 

Kk L = J f^ + ^ ) , J = ±±,±, (67) 

with 

TffiT = |(1 ± V241) , 7^ = ^(-17 ± V241) . (68) 

The non-logarithmic corrections in the basis in which 7g°L is diagonal and LO Wilson 
coefficients Cg LL (j = ±±, ±) are normalized to 1 read 

r s ± L ± L = l(25±v / 241), r± L = \ (7 ± V^41) - (69) 

In the Appendix |A] we list the Wilson coefficients and operators in this new basis. 

In Fig. [6] we show the results for the SLL sector, again for the three different cases: LO 
(dashed blue), NLO with only logarithmic corrections (dotted green) and NLO with both 
logarithmic and non- logarithmic corrections (solid red). As expected from the values of 
i^SLL the dashed blue lines show a very strong /i in dependence in the case of ++ and +. 
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The inclusion of the logarithmic 0(a s ) terms calculated by us practically removes this 
dependence (dotted green lines). However, the large size of the non-logarithmic terms 
r su! + implies significant left-over /i in dependence in these two cases at the NLO level 
(solid red lines) that can only be removed by NNLO corrections. 

SLR case 

At last we discuss the special case SLR where no logarithms appear in the matching 
conditions. Diagonalizing (^/(°) SLR ) T we get 

/nSLR _ nSLR /vnSLR. ^ySLR _ 1 ✓-(SLR _ ° (7n\ 

V+ — V2 " JV Vl , L '+ _ _ jy 1 ' ^SLR — > 

r>SLR n SLR n SLR 1 n SLR , ^ySLR (0)- a at , ^ fn-i \ 

Q- = Q2 , C_ = —C 1 + G 2 , 7 SLR = -67V + — . (71) 

We can read off that the coefficient C+ LR cannot be unity in leading order. Only (17^ LR 
can be normalized to 1 and its [x- in dependence related to the anomalous dimension of 
qSlr j g exac tly cancelled by the /j,- m dependece of quark masses. This means simply that 

7slr = ~2>lm as seen explicitely in Eq. ( 17TT) . 



9 Summary 

In various extensions of the Standard Model (SM) tree level non-leptonic decays of hadrons 
receive contributions from new heavy gauge bosons and scalars. Prominent examples 
are the right-handed W bosons in left-right symmetric models and charged Higgs (-ff ± ) 
particles in models with extended scalar sector like two Higgs doublet models and super- 
symmetric models. Of particular interest are the decays in which the contributing local 
operators involve four different flavours so that penguin operators cannot contribute and 
the decays are simpler to analyse theoretically. Anticipating the important role of such 
decays at the LHCb, KEKB and Super-B in Rome and having in mind future improved 
lattice calculations combined with the QCD factorization approach, we have completed 
the existing NLO QCD calculations of these processes by calculating 0(a s ) corrections 
to matching conditions for the Wilson coefficients of all contributing operators in the 
NDR-MS scheme in any extension of the SM. The main results of our paper can be found 
in Section [51 

Our calculation allowed to reduce certain unphysical scale and renormalization scheme 
dependences in the existing NLO calculations. Our results can also be applied to models 
with tree-level heavy neutral gauge boson and scalar exchanges in AF = 1 decays and 
constitute an important part of NLO analyses of those non-leptonic decays to which also 
penguin operators contribute. 

For completeness we have collected all the relevant formulae necessary to perform the full 
NLO renormalization group analysis that would require the evaluation of the hadronic 
matrix elements in the the NDR-MS scheme. They can be found in Section |6] with all 
ingredients given in Sections |2] and [5j 
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We are aware of the fact that present hadronic uncertainties in the decays considered are 
significantly larger than corrections calculated by us. However, one should recall that 
twenty years ago similar comments were made in connection with NLO QCD corrections 
within the SM [5]. During the last decade after significant progress in lattice calculation 
has been made, the results of the 1990 's constitute an important part of any phenomeno- 
logical analysis of weak decays in the SM. We expect that this will be the case of our 
results in due time. 
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A Change of basis: SLL case 

For our numerics in Section [8] we had to diagonalize the anomalous dimension matrix 
in the SLL sector. In this case the situation is a bit more intricate as four operators 
are involved. The eigenvalues 7sll U = = ' = ) are gi ven m Eq. ( 1681) for N = 3. For 
completeness we list here the Wilson coefficients and operators in the basis in which 
^(o)sll - g diagonal and where LO Wilson coefficients Cg LL (J = ±±, ±) are normalized to 
1. We recall, as found in Section [5j that in the original basis only the coefficient Cf LL is 
non- vanishing in the LO. 



8 SvW^ 1 U^vW^ 2 U 32^41 
± i^=Qa LL , (72a) 



C ±± =2 ( 15 ± >/241)C? LL + C% LL + 2(19 ± \/241)Cf LL + 4(16 ± V241)C% LL , 



- =F -^-] Q SLL + (I * Q SLL + (-L ± _?L_ 

8 T 8^4T;^ + V4 T v / 24T; W2 V 32 ± 32v / 24l 

=F t^=<54 LL , (72b) 

C-sll = J_(2i ± V24T)Cf LL + C 2 SLL + -(1 ± V241)Cl LL + -(-16 =f V2U)Cf LL . 
10 5 5 
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